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A duality transformation that interrelates expanding and contracting cosmological models is 
shown to single out a duality invariant, interacting two-component description of any irrotational, 
geodesic and shearfree cosmic medium with vanishing three curvature scalar. We apply this feature 
to a system of matter and radiation, to a mixture of dark matter and dark energy, to minimal and 
conformal scalar fields, and to an enlarged Chaplygin gas model of the cosmic substratum. We 
extend the concept of duality transformations to cosmological perturbations and demonstrate the 
invariance of adiabatic pressure perturbations under these transformations. 
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I. INTRODUCTION 

Different kinds of symmetries are crucial to obtain so- 
lutions of Einstein's gravitational field equations. One of 
the best known examples are the symmetries that underly 
the cosmological principle which is at the heart of most 
investigations in cosmology. Only the reduction of the 
number of degrees of freedom of the full theory by this 
symmetry principle makes the complicated set of non- 
linear equations tractable and allows one to characterize 
the cosmological dynamics within Friedmann-Lemaitre- 
Robertson- Walker (FLRW) models by a few key param- 
eters, notably the Hubble parameter which describes the 
expansion rate of the observable Universe. The Hub- 
ble parameter enters Friedmann's equation quadratically. 
This gives rise to another symmetry within a class of 
FLRW models itself which will be the focus of this paper. 
Starting point is the observation that, because of this 
quadratic dependence, Friedmann's equation remains in- 
variant under a transformation H —H for the spa- 
tially flat case. This means, it describes both expanding 
and contracting solutions. Such kind of duality between 
expanding and contracting models, except being of gen- 
eral theoretical interest, may be relevant, e.g., to relate 
post-big bang to pre-big bang background and pertur- 
bative solutions in corresponding string theoretical sce- 
narios [H, 0, i, II i, H 0, i, i, The transforma- 
tion H ^ —H can be seen as a consequence of a change 
a ^ 1/a of the scale factor a of the Robertson- Walker 
metric which is called scale factor duality [III [13, [H, [I3l ■ 
This duality becomes a true symmetry of Einstein's equa- 
tions when it is supplemented by the transformations 
yO — !■ p of the energy density and p —2p — p of the pres- 
sure, where we assumed a perfect fluid type structure of 
the stress energy tensor. Under this set of transforma- 
tions Friedmann's equation and the conservation equa- 
tions remain invariant. Interestingly, this internal sym- 
metry comprises both the geometrical quantity H and 
the fluid source quantities p and p. So that, when we re- 
fer to a duality transformation in the following, we have 
in mind the set of transformations H —H, p ^ p and 
p -2p-_p. 



The aim of this paper is to discuss implications for 
the cosmological dynamics under the assumption of an 
invariance of the basic equations under duality transfor- 
mations. We will uncover an underlying invariant two- 
component structure of any cosmic medium with an en- 
ergy momentum tensor of the perfect fluid type. These 
two components have fixed equations of state and are 
converted into each other by a duality transformation. 
One of them is radiation (barotropic index 4/3), the other 
one is of the phantom type and called "dual radiation" 
because its barotropic index is —4/3. These components 
represent reference fluids for the cosmic substratum. Any 
medium with arbitrary equation of state can be decom- 
posed into an interacting mixture of radiation and dual 
radiation. The circumstance that a symmetry require- 
ment, here the duality invariance, is accompanied by the 
appearance of interactions, is a feature that is well known 
from gauge theories. While a decomposition of the cos- 
mic medium into radiation and dual radiation may ap- 
pear artificial, its invariance properties will turn out to be 
useful to study the behavior of physically relevant two- 
component systems like a mixture of matter and radia- 
tion or of dark matter and dark energy. We shall also 
discuss realizations of the duality invariant substructure 
for a constant equation of state, for minimally and con- 
formally coupled scalar fields and for an enlarged Chap- 
lygin gas. For the case that the energy densities of both 
components are equal, we recover the de Sitter universe. 
Finally, we investigate the behavior of cosmological per- 
turbations and show that the concept of adiabatic per- 
turbations is duality invariant. 

There is still another aspect which we would like to 
mention. Present observational data give strong support 
to an almost flat universe. This includes the possibility 
that the universe is indeed exactly flat (fc = 0). Although 
this is an "unlikely" singular case, it deserves attention. 
One may speculate that a so far unknown symmetry en- 
forces the universe to be spatially flat. If e.g., duality in- 
variance of the kind dealt with here is a symmetry which 
is realized in the actual Universe, there would be no flat- 
ness problem. 

The paper is organized as follows. In section |TT] we 
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establish our concept of duality transformations and du- 
ality invariance. In section IIIII we demonstrate that for 
any perfect fluid type cosmic medium the behavior of the 
relevant quantities under duality transformations natu- 
rally singles out a decomposition into a radiative fluid 
and a phantom fluid which we call dual radiation. We 
clarify the relation of this decomposition to different two- 
component descriptions of the cosmic medium. The gen- 
eral formalism is illustrated in section [TVl for a system of 
matter and radiation, for a mixture of dark matter and 
dark energy, for several scalar field configurations and 
for the enlarged Chaplygin gas. Section |V] is devoted to 
the behavior of cosmological perturbations under duality 
transformations. In section lVll we summarize our results. 
Units have been chosen such that Sir G = c = h = 1. 



II. COSMOLOGICAL FLUID DYNAMICS 

Let us consider a cosmic medium with a perfect fluid 
type energy momentum tensor 



the sign of the expansion scalar 



nmn 'tn^,n 



pu u + ph , 



(2.1) 



where p is the energy density, p is the pressure, m™ 
is the four velocity normalized to u"^u„i = —1, and 
j^mn _ gmn _j_ yTa^n ^ y-g^ Einstciu's field equations 



(2.2) 



one has 



together with 



e = -e 



p^p 



-2p-p. 



(2.7) 



(2.8) 



Under this transformation the energy balance (|2.5p re- 
mains invariant due to p + p —t —(p + p). Consequently, 
if the weak energy condition is satisfied in a given cos- 
mological model, i.e. > 0, it is violated in its dual 
and vice versa. In a FLRW cosmology the transformation 
law (|2.7|) implies the transformation rule a — > a = 1/a for 
the scale factor a. Accordingly, if a certain configuration 
(say, the unbarred one) describes a phase of contraction, 
the barred one describes a phase of expansion. These cos- 
mological solutions are said to be dual to each other. In 
particular, there is a duality between a final contracting 
big crunch and a final expanding big rip. In general, the 
equation of state parameter will change under a duality 
transformation according to 



7 



■P 



P 



P + P 
P 



-1 



(2.9) 



where 7 and 7 are the barotropic indices of the fluid 
and its associated "dual fluid", respectively. The only 
invariant case is p = —p. This is related to the fact that 
the de Sitter universe is free of singularities. Using the 
transformation properties (|2.7p for G and (|2.8p for p and 
p in the Raychaudhuri equation (j2.6[) . the latter turns 
out to be duality invariant if 



G"^"'UmUn — P , 



1 



■G""ft„,„=p, (2.3) 



1 



-,Q^^P, 



(2.10) 



as well as 



= 2 (P + 3p) , R = p~3p. (2.4) 



The energy balance reads 

p + e{p + p) = 0, 



(2.5) 



where = u"^^ is the fluid expansion and p = /0,„u". 
For an irrotational, geodesic and shearfree motion the 
Raychaudhuri equation for Q reduces to 



i.e., the three curvature scalar TZ of the hypersurfaces 
orthogonal to u", in the present case given by 
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(2.11) 



has to vanish. Obviously, equation (|2.1ip itself is 
invariant under a duality transformation. The duality 
invariance implies the relations 



= G" 



(2.12) 



e + ie2 + i(p + 3p) = 0. (2.6) 
We introduce a duality transformation by a change in ^ ~ fj^mf^n ^ ^ '^mn^ i^z.io; 
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between the components of the Einstein tensor. Since 
energy momentum conservation is preserved under dual 
transformations, we have also '^iG], = ViG].. Here we 
have introduced the symbol V which is the usual covari- 
ant derivative with the affine connection V\i substituted 
by the corresponding duality transformed quantity F^^ 
(see below). Hence, if the Einstein equations are satis- 
fied for a certain configuration, then a set of equations 
with the same form, in general, describes another poten- 
tially possible cosmological model. 

In a homogeneous and isotropic background the fluid 
expansion in terms of the scale factor a is = iH — 
3a/ a. The equations (|2.6p and (|2.10p reduce to 

= p, and -2H = p + p, (2.14) 

respectively and they are form invariant under the du- 
ality transformation (|2.7p - (|2.8p . In this case the former 
equation takes the form 

H^-H d^a-K (2.15) 

The transformation rule of the deceleration parameter, 
q EE -ti/{aH^) = {1 + 3p/p) /2, reads 

q = -q-2. (2.16) 

Examples are: (i) p = p/3 q = 1 ^ = —3, 
(ii) p = -7p/3 ^ g = -3 g 1, (iii) 

p — 4^ q = 1/2 ^ q — —5/2. In all these cases 
the sign of the deceleration parameter changes, i.e., a 
universe with decelerated expansion is transformed to 
a universe with accelerated contraction and vice versa. 
But for equations of state p < — p/3 both q and q are 
negative, e.g., (iv) p — —2p/3 q ^ —1/2 ^ q ^ —3/2. 
Consistent with the general setting discussed so far we 
have (v) p ~ —p <S4>g = — l=>g = — las the only case 
with q ^ q. The cases (i) and (ii) will be of particular 
importance for the decomposition to be introduced 
in the following section. Finally, for a homogeneous 
and isotropic background, the non-vanishing Christoffcl 
symbols 

F^, = HS^ , r% = a'H5^, (2.17) 

behave as 

-Tqi/ — — Toi/ I f J^j, — j^^uu (2-18) 

under a duality transformation. The transformations 
([2?T2)) and ((27T3)) for the Einstein tensor specify to G^ = 
Gg and Gl^G^ = Gl^ 2G^ - G\. This shows explic- 
itly the invariance of the conservation laws under a dual 
transformation. 



III. TWO-COMPONENT PICTURE, GENERAL 
RELATIONS 

In many periods of the evolution of the Universe the 
cosmic medium is adequately described as a mixture of 
two components. The total energy and the total pressure 
are then split according to 

P^Pa+Pb, P = Pa+Pb- (3.1) 

A typical example is a system of matter and radiation 
which underlies, e.g., simple models of reheating in the 
early universe [15l . Ilq . Another example of current inter- 
est is the description of the present phase of accelerated 
expansion in terms of dark energy and dark matter as the 
two dynamically dominating components of the cosmic 
substratum. Now, the duality transformations (|2.7[) and 
(|2.8p are independent of any underlying two- or multi- 
component description. The behavior of the individual 
components under duality transformations remains open. 
There is a degeneracy with respect to a potential internal 
structure of the medium. What matters there are only 
the total energy density and the total pressure. What 
one would like to have is a prescription of how a given 
two-component structure behaves under a duality trans- 
formation, i.e., which are the counterparts pA and of 
Pa and ps, respectively. It is this question that we are 
interested in in the following. 

A convenient starting point is the circumstance that 
the transformation rule of the scalar curvature (cf. (|2.4p 
and dMl)), 

p-3p ^ p-3p^7p + 3p , (3.2) 

apparently singles out the duality related equations of 
state p — p/3 and p = —7p/3, the latter being of the 
phantom type. Namely, for p — p/3 the curvature scalar 
is zero, while the dual curvature scalar vanishes for 
p = p/3, equivalent to p = —7 p/3. We shall demonstrate 
now that these two (fixed) equations of state, which we 
will denote by 

Pi = y and pa = ~^P2 , (3.3) 

respectively, may be used as a two dimensional basis with 
respect to which any cosmic medium with arbitrary equa- 
tion of state may formally be decomposed. The result is 
an interacting two component system with p = pi + P2 
and p = pi + P2- For this singled out split it will be pos- 
sible to obtain the transformation properties of the indi- 
vidual components 1 and 2, i.e., to remove the mentioned 
degeneracy. Then, in a subsequent step we shall derive 
the dual counterparts pA and pB for any split p = pa+Pb 
of interest, in particular for a system of matter and ra- 
diation as well as a mixture of dark matter and dark 
energy. 
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The curvature scalar and its dual in (|3.2p can be com- 
bined into 



-{R + R) = p = Pi + P2 



(3.4) 



and 



'R-R) = {p + p) = liPi+l2P2, (3.5) 



where 

Pi = (71 - 1) Pi and p2 = (72 - 1) P2, (3.6) 
with 71 = 4/3 and 72 — —71 = —4/3, respectively, and 



R 1 



i? 1 



Pi 



= -(7p + 3p), p2 = - = -(p-3p), (3.7) 



8 8 



8 8 



or 



37 



37 



(3.8) 



P2 = (72 — 1)P2 where 71 = 4/3 = —72- ^ phantom com- 
ponent naturally appears in this split. The specific fea- 
tures of the overall equation of state p = p{p) are entirely 
encoded in the interaction between both components. 
Rewriting Einstein's equations (|2.6p and (|2.10p . we get 



«2 



and 



e 



2 • 



= Pl+ P2 



- 2© = 71 (pi ~ P2) , 



(3.13) 



(3.14) 



respectively. From Eqs. 
it follows that 

Pi = P2, 
P2 



3.6p and (|3.7p together with 



P2 



P2 



-- Pi. 
~7pi. 



where =_pwithpi = (7i-l)pi andp2 = (72- 

While the individual barotropic indices 



(3.15) 
(3.16) 

1)P2. 



with 7 = 1 + p/p (cf. Eq. ([M])). While the first compo- 
nent in (13. 8p is a radiation component, the second (phan- 
tom type) component with a barotropic index 72 = —4/3 
will be referred to as "dual radiation" . It is seen by in- 
spection that an equation of state p — —p is the only 
case in which R — R. Under this condition 11 = —H = 
is valid, i.e., there is a duality between de Sitter and 
anti-de Sitter cosmologies. Obviously, this case implies 
Pi = P2 = p/2, i.e., for this configuration the energy is 
equally distributed on both components. 

Since 7 is related to 71 and 72 via 7p = 7ipi -I- 72P2, 
the balances for the components are 



pi + 67^1 



■IP 



(3.9) 



71 = 71 = 3, 



72 = 72 = -g, 



(3.17) 



and 



remain invariant under the dual transformation (which 
is different from the case investigated in [l3|)j the trans- 
formation behavior of 7 is given by Eq. (|2.9p . 

The relations p.l5p and (|3.16p represent the desired 
transformation behavior for the singled out split p.6p - 
(13. 8p . which will be applied below to physically relevant 
decompositions of the cosmic medium. 

The two-fluid picture of radiation and dual radiation 
components can be seen as a symmetric representation 
of the Einstein equations (|3.13p - p.l4p (for vanishing 
three curvature) under the transformations pi — > p2 
and — > — 0. Local energy conservation under these 
circumstances takes the form 



or 



and 



/92 + 67P2 = -■tT'P I (3.10) 



Pi + O71P1 = ^e^i^ ^7P , (3.11) 
op o 



P2 + 672^2 --^e^-^7P- (3.12) 



P = -3© {Pi - P2) 



and the balances (13. lip and ()3.12p become 



Pi + gQpi = rp2 , 



(3.18) 



P2 - 36/92 = -rp2 , (3.19) 



where 



These relations prove our statement that under the 
condition of vanishing three curvature any irrotational, 
geodesic and shearfree cosmic substratum with effective 
scalar pressure p and energy density p can always be 
described as an interacting two-component mixture p = 
Pi + P2 OLTid p = pi -\- p2 with pi = (71 — l)pi and 



Tp2 = 8H 



P1P2 , 5 f P 



P ■ 



(3.20) 



With the given preferred split into radiation and dual ra- 
diation, the coupling term is completely fixed by the un- 
derlying symmetry requirement. As already mentioned, 
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this reminds of properties that characterize gauge theo- 
ries. However, so far the relation of the expression (|3.20p 
to physically relevant interactions remains open. 
Since the expression for the energy density ratio 



and 



P2 ^ 1 - 3p/p ^ 1 - 37/4 
Pi 7 + 3p/p 1 + 87/4 ' 



(3.21) 



is positive, the total equation of state parameter p/ p is 
consistently found to be in the range —7/3 < p/p < 1/3. 
Obviously, e — is valid. Generally, the dynamics of 
the ratio e is governed by the equation 



(3.22) 



A simple example is a constant equation of state in 
a homogeneous and isotropic FLRW space time. For 
(p/p) = 0, equivalent to 7 = 4(1 — e)/3(l-|-e) = constant, 
equation (|3.22p yields e = const. The energy densities 
scale as 

pi, p2, p (xa~^~ . (3.23) 

The interaction corresponds to a decay of component 2 
into component 1, where the decay rate F, according to 
Eq. |3?20| . is given by 



F = 8if^ 



8H- 



1 



(3.24) 



i.e., the rate F is proportional to the Hubble rate. Al- 
though the component 2 decays, the ratio e = P2/P1 re- 
mains constant. This is due to the circumstance that 
72 is negative, i.e., without interaction p2 would grow 
with the expansion (for H > 0). Growth and decay are 
balanced such that the ratio of the energy densities is 
unaltered. 

Interesting subcases are: (i) radiation, 7 = 4/3, equiv- 
alent to pi = p and e = 0, (ii) dust 7=1, where 
pi = 7p/8 and e = 1/7, (iii) vacuum energy, 7 = 0, 
characterized by pi = p/2 and e = 1, (iv) dual radiation, 
7 = —4/3, for which pi = and e — > 00 . 

Now, there will be hardly an epoch in which the decom- 
position into radiation and dual radiation with a specific 
interaction between both components is a physically ad- 
equate description of the cosmic medium, although it is 
formally always valid. However, the duality invariance of 
this split can be used to characterize the duality behav- 
ior of other, more realistic two-component descriptions, 
which by themselves are not duality invariant. 

Let us consider to this purpose a different decomposi- 
tion of the cosmic medium according to (|3.ip with 



(3.25) 



PA = ilA - 1) PA, PB = (7-B - 1) PB 



The relations between the duality invariant split p. 71 
and the arbitrary decomposition (|3.ip are 



Pi 



[{■i + 3ja) PA + i-i + '-^Ib) Pb] 



(3.26) 



P2 = o [(4 - 37a) pa + (4 - 375) pb] 



and 



PA 



PB 



7B - I 

IB - lA 



Pi 



7B + I 
7s - lA 



P-2 



7A-| 



_1A ~ IB 

Of course, this satisfies 



Pi 



lA 



lA - IB 



■P2 



Pi + P2 = PA + PB 



(3.27) 



(3.28) 



(3.29) 



(3.30) 



while the difference pi— p2 which appears in (|3.18p . trans- 
forms into 



Pi- P2 = ^ [lAPA + 7i3PB] 



(3.31) 



The transformation properties (|3.15p and p.l6p of pi , p2 
and pi, p2, respectively, can now be used to obtain the 
corresponding behavior of pA, ps, PA and pB- In the 
following we discuss several applications of this scheme. 



IV. TWO-COMPONENT PICTURE, 
APPLICATIONS 



A. Matter and radiation 

As a first application we consider a mixture of matter 
(subscript to) and radiation (subscript r). 



n 

= , P-'^ ^ 



(4.1) 



With 



7A = 7m = 1, 75 = 7r = - (4.2) 



we obtain from ([3:281) and 

Pm = 8p2, Pr^ Pi- 7p2 , (4.3) 

or, from (|3:26| and ([3^, 

Pi = -^Pm +Pr, P2 = ■ (4.4) 



It follows that 



Pi - P2 = ^Pm + Pr • 



(4.5) 



The behavior under duality transformations is 

Pm = 7pm + 8pr, pr = -6pm - 7 Pr , (4.6) 
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guaranteein g + Pr = Pm + pr- 
Equation p.lSp can be written as 

Pm + iHpm + Pr + -iH Pr = , 

which is equivalent to 

p,n + ZHprn = —Q , 



and 



where 



Pr + 4:Hpr = Q , 



g = (r - 7H) p„ 



(4.7) 

(4.8) 
(4.9) 

(4.10) 



The set of equations (|4.8p - (|4.10p describes matter in 
interaction with radiation. While T characterizes the in- 
teraction between radiation and dual radiation and may 
be considered a fictitious interaction rate, the quantity Q 
is supposed to describe a real physical interaction, e.g., 
the decay of a scalar field into radiation in simplified 
reheating models of inflationary scenarios [isl . [l^ . How- 
ever, only the knowledge of the duality transformation 
behavior —Tp2 of p.20p allows us to obtain the 

dual counterpart Q oi Q. The combinations 



and 



Q = 56H (pi - p2) 



-Q = 16rp2-56H{pi+p2) 



(4.11) 



(4.12) 



are obviously duality invariant. The quantity Q itself, 
however, changes under a duality transformation. From 
(|4.10p we find that the non-interacting case Q = cor- 
responds to r 7H. Its dual Q, on the other hand, 
is different from zero. Consequently, the dual of a non- 
interacting mixture of matter and radiation is interact- 
ing. From 



1_ 
-fPr 



~2p„ 



and 



Pr, 



Pr 



7 



Pr) 



(4.13) 



(4.14) 



we obtain the corresponding transformed equations of 
state 

1 



P,n =0, Pr = -pr , (4-15) 

i.e., the individual equations of state preserve their form. 



B. Dark matter and dark energy 

Another interesting special case is a mixture of matter, 
7m = 7^ = 1 and dark energy 7x = 7_b = 0. Then (cf. 



pi 



;Px 



7 



-yPM 



(4.16) 



and 



P2 = -j^Px + -^PM 



PM = 3 (/'I - P2) 



Px = ~3^i + 3/^2 



(4.17) 



(4.18) 



(4.19) 



With these transformations the balance equation p.l8p 
takes the form 



Pm + iHpM = -px 



(4.20) 



This represents a model in which matter and dark energy 
are in general interacting with each other: 



pM + iHpM = Q, 



PX = 



(4.21) 



with 



Q - -^H Upx + 1pm j + (px + \pM ] ■ (4.22) 



Under duality transformations the variables of this model 
behave as 



PM 



The combinations 



-PM, Px = Px + "^PM 



Q^8H (pi - P2) 



and 



y [Tp2-2H{p,+p2))] 



(4.23) 



(4.24) 



(4.25) 



are duality invariant. 

The ACDM model is realized for Q = 0. Its dual 
counterpart is an interacting model with 



Q^8H (P1-P2) 



H pi+ 7p2 

2 P2 



(4.26) 



Any case Q ^ corresponds to an interaction between 
dark matter and dark energy. Models of this type have 
attracted considerable attention in the literature, in par- 
ticular with respect to the coincidence problem '18', 13| 

Similarly to the previous example, the transformed 
equations of state are obtained via 



The result is 



Pm + Px 



PM = 0, 



{P- 



'M + Px ) 



Px 



'PX ■ 



(4.27) 



(4.28) 



Again, the equations of state are form-invariant while 
the interaction rate in the original model differs from the 
corresponding rate in the transformed model. 
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C. Systems with one constant component 



The minimal scalar field 



The balance equations (|3.18[) - (|3.20[) simplify if one 
of the components remains constant. For p2 = 0, e.g., 
it follows from (|3.19p that T — AH . This corresponds to 
the case 7a =4/3 and 73 in (|3.31|) . Then pi — p2 = 
Pa = N (X. a^'^ and pg = K = const. Similarly, for 
pi = one has F = 4H pi/p2 and p2 — Pi = F cc a'^. 
For these configurations the radiation and dual radiation 
components can be written as 



pi = iV- 



A 
2"' 



A 



P2 = 



or 



Pi = TT' 



where p = N 



A 

A or p 



A 



(4.29) 



(4.30) 



F + A and A is constant for 
the given dynamics. However, A changes under a duality 
transformation. The quantities N{F) and A transform 



N{F):^~N{F), A^2N{F)+A 



(4.31) 



Obviously, Eq. (|4.29p describes a universe filled with radi- 
ation and a cosmological constant, while Eq. (|4.30p is an 
analogous configuration with radiation replaced by dual 
radiation. 

Because of the mentioned dependence on the scale factor 
N and F can be written as 



N 



No 



F = Fo a\ 



(4.32) 



where iVo and _Fo are integration constants. In terms of 
iV (or F) and A the Friedmann equation (|3.13p becomes 



No 



A, or 3i?2 = i^o a* + A 



(4.33) 



Integration of equations (|4.33p provides us with the scale 
factors 



''N 



±,/^sinh^^t 



^iVo , /4A^ 
— — cosh \ / — t 
A V 3 




iVo > , (4.34) 



iVo < , (4.35) 



in the radiation case and 



a^iNo 



Fo 



ap 



aN{No = Fo) 



(4.36) 



for the dual radiation case. The solution ap, defined in 
the region t < 0, describes a big rip at t = 0. On the other 
hand, the super-accelerated solution ajv has an infinite 
time span and a final de Sitter stage. 



In this part we derive an explicit scalar field dynam- 
ics which is a realization of the structure (|4.29p - (|4.30p 
with the results (I4.34p - (|4.36p . Starting point are the 
expressions 



1 



1 



r + V{(f>) and p^^ -,p' -Vi(t>) (4.37) 



for the energy density p^ and the pressure p^, respec- 
tively, of a scalar field. The decomposition into radiation 
and dual radiation according to Eq. (13. 7p then yields 



1 



Pi 



92 



(4.38) 



where the first component represents some kind of 
quintessential scalar field while the second component 
is associated with a phantom scalar field. Independent 
of any interaction between the two fluid components, the 
energy-momentum tensor conservation of the system as a 
whole is equivalent to the Klein-Gordon equation. Apply- 
ing the transformation rules (|2.8p to p0 and p,/, of (|4.37ll . 
we obtain the transformation properties for the kinetic 
energy and for the potential of the scalar field: 



i2 



V - 



12 



(4.39) 



To see how the structure ()4.29p - ()4.30p applies to the 
scalar field dynamics, it is convenient to choose the po- 
tentials 



Vs^^(j)'^+A and 



5 



A, 



(4.40) 



in (|4.38p . As will become clear in a moment, Vs and Vph 
are associated with "conventional" scalar (subscript s) 
and phantom (subscript ph) fields, respectively. Namely, 
with (|4.40p the components pi and p2 in (|4.38p take the 
form 



Pis 



A 
2 



and 



P2s 



A 

2' 



respectively, or 
A 



Piph = and 



3 



P2ph 



A 
2" 



(4.41) 



(4.42) 



These components have indeed the structure of (|4.29p 
and (|4:30l) with N = F^ = 30V4 = -Fph = -F. There- 
fore, the corresponding scale factors are given by (|4.34p - 
(|4.36p . Inserting the non-phantom and phantom fluids 
(|0T|) - (|0^ into the balance equations ([?TT1) - (|XT^ and 
integrating, we obtain the first integrals 



4A^o _2 



and 



Pph 



4^^o 2 

a 

3 



(4.43) 
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for the scalar (non-phantom) and phantom field cases, 
respectively. The explicit scalar field dynamics is then 
derived by using Eqs. ([434| - ([436| in (j443l) . For the 
potentials (|4.40p we find 



K = A 



cosh 



■ sinh 



with 



and 



with 



or, 



Intanh^/ — Nq > 0, 



= A 



cos'^ ~\ — sin^ iq 
3 



(4.44) 



(4.45) 



(4.46) 



= -2iarctaneV4A73t^ TVq < 0, (4.47) 



ph 



A 



cosh i6 + 



■ sinh^ ia 



with 



and 



with 



= — i Intanh ^/ — Fq > 0, 



cos^ d) sin^ i 



: 2arctane\/^', iVo < 



(4.48) 



(4.49) 



(4.50) 



(4.51) 



We mention again, that the duality related components 
pi and p2 in Eq. (|4.38p describe a quintessence-type 
scalar field and a phantom field, respectively. This 
is quite similar to features which are characteristic for 
"quintom" models of the cosmic substratum . A mix- 
ture of this type has attracted considerable interest re- 
cently since it may describe the crossing of the so called 
phantom divide which seems to be favored by a number 
of observational data l2l|, [S^l . 



2. The conformal scalar field 

An interesting particular case is obtained for a van- 
ishing N, corresponding to the one-component limit 
Pi — P2 — PA = 0- Namely, for = the constant A 
preserves its form, i.e., A = A (cf. Eq. ()4.31|) ). We pro- 
vide a representation of this case by a conformal scalar 
field i/j with the potential V{ip) = Xip'^ + A with A > 0. 



This potential has received much attention in the litera- 
ture in connection with the early inflationary epoch [23l | . 
This simplified model leads to a final accelerated expan- 
sion phase retaining the essentials of minimally coupled 
approaches. The energy density and the pressure of the 
conformal scalar field are 



P^-^(^ + ^V')' + AV/ + A, 



A. 



(4.52) 
(4.53) 



The radiation dual radiation decomposition p.7p of this 
conformal scalar field leads to 



Pi 



1 r 
2 



(?/> -I- HiPf + 2Xil;^ 



A 

2' 



P2 



-A. (4.54) 



Comparing the first equations in (|4.29p and (|4.32p with 
(|4.54p . we get 



i(^ + iJ^)2-FAV/ = 5, 



(4.55) 



and its corresponding dual (barred) counterpart. Ap- 
plying the transformation rule p.l5p to the components 
(I4.54P , we obtain 



+ H^jf + W = -^(V^ + H^f - , (4.56) 
A = A + (?/> + Hi^f + 2AV'''. (4.57) 



Now we compare Eqs. (j4.55p . its dual (barred), and 
It follows that No/a^ = -Nq/o^ ^ N^a^. This 
has a solution Nq ~ Nq — 0, for which A = A (cf. 
Eqs. (i43T|) and ([432]) ). Consequently, iV = -iV = 0, 
which is just a realization of the special situation men- 
tioned above. In this case the Friedmann equation sim- 
plifies to ?>H^ = A, w ith the de Sitter solution a — e^"*, 
where Hq = ±a/A/3. 

Integrating Eq. (|4.56p we obtain the conformal field 



V^[i 



A < 0, 



(4.58) 



and its transformation rule 



Hq 



(4.59) 



So, the duality invariance of the de Sitter solution induces 
a linear transformation group acting on ij:. In the par- 
ticular case of A = A, i.e., for an invariance of the field 
strength, the radiation - dual radiation decomposition 
produces a simple translation of the conformal field. 

This completes our discussion of the internal structure 
generated by our duality invariant two-component de- 
scription when applied to a (minimal or conformal) scalar 
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field dynamics. It is expedient to notice that the under- 
lying decompositions (|4.29p and (|4.30p are relations in 
which the energy densities of the components, pi and p2, 
depend linearly on the total energy density p. In the 
following subsection we shall perform a fluid dynamical 
analysis for an enlarged Chaplygin gas, for which the cor- 
responding decomposition is non-linear. 



D. Enlarged Chaplygin gas 

So far we have emphasized that a given cosmological 
model can be decomposed with respect to a duality in- 
variant two-component basis. Here we focus on the in- 
verse procedure. Assuming a fairly general non-linear 
structure for pi and p2 as a starting point, we shall 
construct an effective equation of state for the cosmic 
medium as a whole. More specifically, we shall derive 
in this way a set of cosmologies (first considered in 241) 
which were called enlarged Chaplygin cosmologies in 25 1 . 
This class includes as particular cases generalized 26|, 



It results in an energy density 



extended and modified Chaplygin cosmologies [l^l along 
with their duals. This larger class of Chaplygin cosmolo- 
gies, which can be induced using duality transformations 
[25| , represents a further generalization of the extended 
and modified Chaplygin gas scenarios, because it includes 
both super-accelerated and contracting cosmologies. 

We start by introducing general radiation and dual ra- 
diation components 

Pi=bp + J{p) and p2 = {l-b)p-.J{p), (4.60) 

respectively, where 6 is a constant and J{p) is an 
arbitrary, in general non-linear, function of the total 
energy density. Of course, p = pi + P2 is valid. The 
transformation rules (|3.15p imply 



J+J+{b + b)p^p 



(4.61) 



P ■ 



A 



B 



a 



37o(1-I-q;) 



1/1+a 



(4.66) 



with an arbitrary integration constant B. It is conve- 
nient to introduce the quantity y = Aa^'"'^^^°'^ / B . The 
duality transformations (j2.8p require that 70 — —70 and 
a — 1/a, which guarantees the invariance of y. This 
behavior is also consistent with the transformation rule 
p.9|) of the overall barotropic index 7 = 70/(1 + y). For 
the time dependence of the energy density ratio e we 
find 
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(1 + 2/r 



(4.67) 



For any H > and y > the ratio e, given by (|3.2ip . 
increases. For large values of the scale factor, correspond- 
ing to large y - values, the quantity 7 tends to zero, i.e., 
p — p, e — > ef = 1 and the scale factor tends to the 
stable (see below) de Sitter solution. 

For small values of the scale factor, equivalent to small 
values of y, i.e., at an early cosmological epoch, we 
have 7 — > 7o and e tends to the constant initial value 
£i — ^ilo)- Remarkably, all these limits are independent 
of a. There is an evolution of e from < 1 to e = 1. 
For the latter value, equivalent to p = — p, we recover 
the expected result that the energy densities of radiation 
and dual radiation coincide. The evolution of the density 
ratio reminds of a similar feature in interacting cosmo- 
logical models in which a suitable interaction is used to 
address the coincidence problem [2^. The decay rate 
(cf. Eq. (|3.20p ) by which energy is transferred from com- 
ponent 2 to component one is 



AH 1 



4' 4^ ' 1 



37/4 



(4.68) 



For this configuration the overall barotropic index 7 — 
4(pi - P2)/3p is 



(4.62) 

while the equation of state of the total fluid is given by 

- = ^(86-7) + !^. (4.63) 
p 3 3p 

Now we choose 



5=1 + ^ and J = -^ (4.64) 
2 8 8p" ^ ^ 

with constants 70, A and a. By this procedure we are 
led to the equation of state for the enlarged Chaplygin 
gas. 



- =70-1-70— Y 
p p"+^ 



(4.65) 



It changes between r„„iy = 4i?(l -t- 37o/4) for y ^ 1 at 
small values of the scale factor to Fi^t^ = 4iJ for ?/ ^ 1 
in the late time limit with e — Ci — 1. These limits of 
the decay rate coincide with those obtained from (|3.20p . 
So, the enlarged Chaplygin gas interpolates between a 
scaling era and a de Sitter scenario. 

To investigate the stability properties of the stationary 
solution we consider e — + 5 with eo = and assume 
that 5 is an arbitrary function. Inserting this expression 
for e in Eq. (|3?2T|l we find 



24 



(4 + 37o)- 



■7. 



(4.69) 



where 



7 = -3ff (1 + 0)70^^^^ . (4.70) 



There are two cases: 
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a) If (5 > with S <0 and 7 > 0, then (1 + a) y < and 
the constant solution eg is stable. 

b) If (5 < with 6 > and 7 < 0, then (1 + a) y > and 
the constant solution eg is stable as well. 

Case b) covers the physically most interesting configura- 
tion A > 0, B > Q, a > Q. In particular, the stability 
property holds for the asymptotic solution eg = £/ = 1- 

Finishing this part, we mention that in the literature 
there exist other types of decompositions of Chaplygin 
gases which describe two-component mixtures of dark 
matter and dark energy [2^, [3d|. These different splits 
(and their behavior under duality transformations) may 
be recovered from the present split into radiation and 
dual radiation in a similar way in which, c.g, the ACDM 
model was obtained in subsection IIVBI 



A. Matter perturbations 

To study perturbations about the homogeneous and 
isotropic background we introduce the following quan- 
tities (cf. [13] )■ We define the total fractional energy 
density perturbation 



D ■: 



P+P 



-3H- 



(5.6) 



and the energy density perturbations for the components 
(A =1,2) 



Da = -3H^ 

PA 



(5.7) 



PERTURBATION DYNAMICS AND 
DUALITY 



(The subscript (&) will be omitted from now on.) Under 
a duality transformation these quantities behave as 



So far we have studied implications of a duality 
motivated split of the cosmological dynamics under 
the conditions of spatial homogeneity and anisotropy. 
However, the formalism outlined in sections HIl and [1111 is 
more general. It allows us to split the entire dynamics 
into a homogeneous and isotropic background part, 
described in section IIV( and first-order perturbations 
about this background. All the quantities p, p, 8 are 
assumed to decompose according to 



P = P(b) + P, P = Pit) + P, = 0(6) + 0, ■• 



(5.1) 



where the subscript (6) stands for background. We re- 
call that the duality transformations (P?7)) . (PTH]) . (|XT51) 
(|3.16p . (|3.17p were supposed to hold for the total quan- 
tities p, p, 0, etc. On the other hand, all the previous 
applications relied on the fact that these relations are 
separately valid in the background, i.e.. 



P(b) - P[b), P{b) = -2p(fc) - p(f,). 



0(6) — ^0(6); 



(5.2) 



Perturbations of the duality transformed quantities are 
introduced by 



P^P(b)+P, P^P(b)+P, 



= 



(b) 



0, 



(5.3) 



Taking here 
tion behavior 



into account leads to the transforma- 



D = -£>, Di = -D2, D2 = -Di 



(5.8) 



The fractional energy density perturbations of compo- 
nent 1 transform into corresponding perturbations of 
component 2 with a reversed sign and vice versa. The 
total energy density perturbations just change their sign. 
The pressure perturbations are 



P 



P 



-3H 



P 



-3H 



pp 
PP 



and 



P^^_3i/^ = -3iI^^. 



PA 



PA PA 



They transform as 



P = P + 2D, Pi = -P2, P2 = 7Pi 



One also realizes that 



and 



(5.9) 



(5.10) 



(5.11) 



(5.12) 



p = p, p = -2p-p, e = -0. 



(5.4) 



P P 



(5.13) 



of the perturbations. Also, the relations 



Pi=P2i P2 = Pii Pi=-yP2, P2 = (5.5) 



are valid. 



While the individual quantities Pa/ Pa Siie duality 
invariant. 



El - il - 1 El - El - -1 
pi Pi 3' p2 h 3 ' 



(5.14) 
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the total adiabatic sound speed is not. It transforms as 

(5.15) 



P 
P 



P 



Again, an equation of state p — — p is the only case of an 
invariant "sound speed" . Combining (|5.6p , (|5.1ip and 
(|5.15p it follows that the notion of adiabatic perturbations 
is duality invariant: 



P = tD 



4^ 



p = Pd 



p p 

Non-adiabatic pressure perturbations 
P 



P P 
are generally characterized by 

P 



P _ P 
P P 



P 



:D = 



P 





_ in 


-) 




pi 










p2 


.P2 


Pi. 



P2 



7^0 



P2 - 



(5.16) 



(5.17) 



P 

[D2 - Di] 



El 

P2 



(5.18) 



The first two terms on the right-hand side describe 
internal non-adiabatic perturbations within the indi- 
vidual components. The last term takes into account 
non-adiabatic perturbations due to the two-component 
nature of the medium. It is straightforward to check 
that also the non-adiabatic pressure perturbations are 
duality invariant, 



P - ^ P ~ -D 



(5.19) 



Since in our case both components are adiabatic on their 
own (recall that pi = pi/3 and p2 = —7/52/3), the first 
two terms on the right-hand side of (|5.18p vanish identi- 
cally. 



B. Metric perturbations 

A homogenous and isotropic background universe with 
scalar metric perturbations and vanishing anisotropic 
pressure can be characterized by the line element (longi- 
tudinal gauge, cf. [sU p 



(1 + 2^) dt^ + (1 - 2^-) 6a/3dx"dx^ . (5.20) 



The perturbation dynamics is most conveniently de- 
scribed in terms of the gauge-invariant variable [32l | 



-iJ- 



1 



2, p + p 



(5.21) 



Corresponding quantities for the components are 

Ca~~^-H^. (5.22) 

PA 

On large perturbation scales the variable C obeys the 
equation (cf [33i, ,34,. i35]) 



C = -H\P--^D 



(5.23) 



The definition (|5.2ip is motivated by the circumstance 
that under infinitesimal coordinate transformations 

x"' =x"-r(x) (5.24) 

the quantities ij) and p behave as 

^'^^-He, p'^p + p^- (5.25) 

Obviously, it makes sense to define dual metric pertur- 
bations by 



■0 = —1p 



c = -c, 



where 



(5.26) 



(5.27) 



For dual adiabatic perturbations the quantity <^ is ap- 
proximately conserved in the same sense in which the 
variable ^ is conserved for usual adiabatic perturbations. 
For the components we have 

Ci = -C2, C2 = -Ci, (5.28) 
with the consequence that 

C2-Ci=C2-Ci. (5.29) 
Under these circumstances Eq. (I5.23P takes the form 



c 



-3H 



PlP2 



El 

P2 



El 
Pi 



[(2 - Cl 



(5.30) 



where 



El - El 
P2 Pi 



(5.31) 



The simplest way to write an equation for 5' = C2 — Ci 
is in terms of the effective interaction pressure 
U= -rp2/{3H) (cf. Eqs. (I3l9l) and JSSQi). Obviously, 
this quantity is invariant under duality transformations, 
i.e., n = n. The equation is 



S = 3HII- 



P1P2 



(5.32) 
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The gauge invariant quantity Cn describes the perturba- 
tion of the interaction term. It is defined in analogy to 
the other perturbation quantities, 



Cn^-^-H-^, n = -^. (5.33) 
From (|2.1ip we find the perturbed three curvature 




(5.34) 



Since on large perturbation scales (upon neglecting 
spatial gradient terms) 



Under the given conditions we have p2 = epi + epi and 
the quantity (2 — Ci on the right hand side of Eq. I5.30p 
becomes 



C2 - Ci = ^ = 




(5.38) 



Given this expression, Eq. I5.32p then fixes the time de- 
pendence of e. 

With (|5.38p equation (|5.30p for ( simplifies to 

C = 2i/-^, (5.39) 

1 — 

or, in terms of {T/H), 



p e 



(5.35) 



holds, we obtain with (|2.10p the (duality invariant) result 
TZ = 0, i.e., vanishing three curvature perturbations. 

Generally, the behavior of perturbations under a trans- 
formation between contracting and expanding phases will 
be of interest if one wants to trace back primeval post- 
big bang perturbations that are responsible for structure 
formation in our Universe to inhomogeneities in a con- 
tracting pre-big bang phase P, i, i, i, i, H, S i, i, [lO] ■ 
While we do not aim to study specific scenarios of such 
a type in this paper, we shall briefly characterize a sim- 
ple toy model for non-adiabatic perturbations and see 
whether and how it exhibits the expected behavior. 

To this purpose we consider the special case e = 
but we allow for spatial variations e ^ (cf. J30J). This 
is a simple way to equip the medium with an internal 
structure. A fluctuation of e ^ corresponds to a 
fluctuation of the equation of state parameter p/ p, 



3 (1 



(5.36) 



Alternatively, using the interaction rate F in (|3.24p . we 
may write 



(5.37) 



which relates the fluctuations of the equation of state 
parameter to fluctuations of the ratio between the inter- 
action rate and the Hubble rate. Fluctuating interaction 
rates are known to produce curvatu re p erturbations in 
certain inflationary scenarios [s^, [s^, Haj . Here, in a dif- 
ferent context, they are used to characterize an internal 
structure within the cosmic substratum ISOl. l39| . 



H /r 



37 \H 



(5.40) 



This demonstrates that e 7^ 0, equivalent to (T/H) ^ 0, 
is necessary to have a non-adiabatic contribution at 
all. Only a fluctuating e can give rise to non-adiabatic 
perturbations. Adiabatic perturbations are character- 
ized by {T/H) = and C = const. For {T/H) = the 
two-component picture does not add any new feature 
to the perturbation dynamics of the single component 
description. With e — e the transformation property of 
the perturbed ratio T/H is 



(5.41) 



The equation (|5.40p preserves its form for H —H, 1 -I- 
p/p ^ -1-p/p, (T/H)^ - {T/H), C ^ -C- Obviously, 
the case 7 = l+ p//9 = 0, equivalent to e = 1, is singular 
and requires e to vanish. Of course, a toy model with 
e = cannot describe a dynamical evolution to e = 1. 
Nevertheless we believe this simple model to be useful to 
demonstrate some general features of the perturbation 
dynamics under duality transformations. 

Wi th (l5?33l) . where H = (72 - 7i) P2/ (1 + e), and 
(|5.38p we obtain 



(5.42) 



for the flrst term in the bracket on the right hand side 
of Eq. I5.32p . Since the second term is 



(5.43) 
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it follows that 

S^O ^ S = const , (5.44) 

i.e., our toy model is characterized by a constant value of 
the isocurvature quantity S. A constant S implies that e 
and {T/H) are constant as well. Notice that S — S (cf. 
Eqs. ([Qg]) ) is consistent with Eqs. ([OS]) and ([OT|) . 

For a constant (T/H) equation (|5.40p may be inte- 
grated to yield 

where the subscript i denotes some initial value. A fluctu- 
ating interaction rate causes a variation in the curvature 
perturbation which is logarithmic in the scale factor. 

VI. DISCUSSION 

Spatially flat cosmological models of General Relativ- 
ity are invariant with respect to a duality transformation 
that converts expanding into contracting scenarios and 
vice versa. We have shown that this property singles out 
a duality invariant decomposition of any cosmic medium 
into two components with fixed equations of state. One 
of them is radiation, the other one is a phantom fluid, 
called dual radiation, characterized by a barotropic in- 
dex —4/3. The specific features of the cosmic medium 
with a generally time depend total equation of state are 
encoded in the interaction between radiation and dual ra- 
diation. The appearance of internal interactions within 
the cosmic medium as the result of a symmetry require- 
ment may be seen as an analog to corresponding prop- 
erties of gauge theories. On the basis of the split into 



radiation and dual radiation we derived the behavior of 
a system of matter and radiation and of a mixture of 
cold dark matter and dark energy under duality trans- 
formations. In general, the duality transformation of a 
non-interacting mixture of two cosmic fluids will result in 
an interacting two-component system. Further examples 
of this decomposition are a fluid with constant equation 
of state, minimally and conformally coupled scalar fields 
and an extended Chaplygin gas. Both scalar field con- 
figurations belong to the subclass for which one of the 
components has constant energy density. The enlarged 
Chaplygin gas represents an example for a non-linear de- 
pendence of the component energy densities on the total 
energy density. The only duality invariant equation of 
state is p = —p, i.e., a de Sitter universe, dual to its anti 
de Sitter counterpart. For this case the energy densities 
of radiation and dual radiation are the same, i.e., each of 
them contributes a fraction of one half to the total energy 
density. 

Finally, we provided an outline for the behavior of 
cosmological perturbations under duality transforma- 
tions. In particular, we demonstrated that the concept 
of adiabatic perturbations is duality invariant. For a 
simple toy model we also calculated a non-adiabatic 
contribution to the curvature perturbation which is char- 
acterized by a logarithmic dependence on the scale factor. 
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